Observation of the decay T(45) — > Y(15)7r + 7r 
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We study transitions between T states with the emission of charged pions using 477 fb^ 1 of data 
collected with the Belle detector at the KEKB asymmetric-energy e + e~ collider. We select inclusive 
T(4S) — > n + fi~n + n~ X events (where X represents anything) and observe a peak in the distribution 
of the mass difference AM = (M M+M -^ +7r - -Af M+M - ). This peak, at AM = (1119.2 ± 0.4) MeV/c 2 , 
is identified as a signal for the decay T(4S) — + T(lS)ir + n~ with a subsequent transition T(IS') — > 
fj, + fi~ . The measured product branching fraction B(T(4S) — > T(IS)iv + iy~) x S(T(15) — ► fi + fi~) = 
(4.4 ± 0.8(stat.) ± 0.6(sys.)) x 10" 6 . When the PDG value for B(T(1S) -> H + fJ,~) is used, this 
corresponds to B(T(4S) -> T{1S)-k + tt~) = (1.8 ± 0.3(stat.) ± 0.2(sys.)) x 10~ 4 and a partial 
decay width r(T(4S) -> T(lS> + 7r-) = (3.7 ± 0.6(stat.) ± 0.7(sys.)) keV. 



PACS numbers: 13.25.Gv, 14.60.Ef, 14.40.Aq 

The bottomonium state T(4S I ) has a mass above the 
threshold for BB pair production and decays mainly 
into 5-meson pairs (B(T(4S) -> BB) > 96% [J). 
However, the decay modes T(AS) — > T{mS)itn with 
m = 1,2 should exist. The first preliminary evidence 
for these decays was presented in Ref. Q, and recently 
the BaBar Collaboration has published an observation of 
these modes Q. A similar decay mode of the charmo- 
nium state ■0(3770) has recently been observed In 
this paper we report the observation of the decay mode 
T(45) -> T(15)7r+7r- from the Belle experiment. The 
results reported here supersede those of Ref. Q • 

We use 477 fb" 1 of data collected on the Y(4S) reso- 
nance and in the nearby continuum to study T(AS) — > 
T(lS')7r + 7r _ decays with a subsequent T(1S) — > fi + fx~ 
transition. Charged particles are reconstructed and iden- 
tified in the Belle detector [||, which consists of a cen- 
tral drift chamber (CDC), aerogel threshold Cherenkov 
counters (ACC), time-of- flight (TOF) scintillation coun- 
ters, an electromagnetic calorimeter (ECL), and a Kl- 
muon detector (KLM). We require that charged tracks 
be well-measured and consistent with originating from 
the interaction point. Charged particles are assigned 
a likelihood Hi [6j (i = fi, n, K) based on the match- 
ing of hits in the KLM to the track extrapolated from 
the CDC, and identified as muons if the likelihood ratio 
P M = Cfj,/ (C^+Ck+Ck) > 0.8, corresponding to a muon 
detection efficiency of approximately 91.5% over the po- 
lar angle range 20° < 8 < 155° and the momentum 
range 0.7GeV/c < p < 3.0GeV/c in the laboratory 
frame. Electron identification uses a similar likelihood 
ratio P e based on CDC, ACC, and ECL information. 
Charged particles that are not identified as muons and 



have a likelihood ratio P e < 0.05 are treated as pions. 

Identification of 7's is based on information from the 
ECL. Calorimeter clusters not associated with recon- 
structed charged tracks and with energies greater than 
50 MeV are considered as 7 candidates. 

Candidates for T(45) — > T(lS')7r + 7r _ decays with the 
subsequent leptonic decay T(15) — * are selected 

from the standard Belle hadronic event sample. The 
most important selection criteria for this event sample 
are the following: multiplicity of charged tracks in an 
event -/V c h > 3; the event's visible energy E vis > 0.2^/s, 
where yfs is the center-of-mass (cm.) energy; the sum 
of good cluster energies in the ECL must satisfy 0.1 < 
^sum/V* — 0.8; the sum of the z components of each 
charged track's and photon's momenta is required to sat- 
isfy \P Z \ < 0.5^/s. The variables E v - ls , E sum , P z are eval- 
uated in the cm. system. 

To select T(15) — * /i + /J~ decays, hadronic events 
are required to contain a pair with Mu+„- > 

9GeV/c 2 and also to satisfy 10.5 GeV < £^ b < 
12.5 GeV, where E!^ is the event's visible energy calcu- 
lated in the laboratory frame. The latter requirement re- 
duces background from poorly reconstructed events. Af- 
ter these requirements, 1.32 x 10 5 events remain. We 
then require the presence of a ir + ir~ pair. To reduce 
background from T(IS') production in radiative return 
processes Q with the subsequent conversion of the emit- 
ted photon into an e + e~ pair that is misidcntified as 
7r + 7r~, we impose an additional requirement on the an- 
gle between the pion momenta in the laboratory system: 
cos^n-Tr < 0.95. The number of selected events is 1084. 

To observe resonance states that decay into the 
T(1S) tt + it~ final state, the distribution of the mass 
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FIG. 1: The mass difference AM = (M„ +/1 - x+7r - -M M+M _) 
distribution where M„+„- hes in the T(1S) mass region. Ar- 
rows show the positions of the mass differences (rn r ^ 2 .s) — 
mrps)), (m T (3S) - m T (is))) and (m T(4S) - m T(ls) ), respec- 
tively, based on PDG values. 




1.2 

AM (GeV/c 2 ) 



FIG. 2: The fit to the third peak in the AM distribution 
(|M M + M - — mx(is)\ < 60MeV/c 2 ) using a Gaussian for the 
signal and a third order polynomial for the background (dot- 
ted line). The solid line shows the sum of a Gaussian and a 
polynomial function. 



' fj,^ [1~ 7T" r 7r - 

m T(lS)l 



Mp+n-) with M t 
< 60MeV/c 2 



is ex- 



diffcrence AM = 
in the range \M^ 

amined (see Fig. [1} for the on-resonance sample. Here 
m r(is) is the nominal T(15) mass. Three peaks are 
seen in the AM distribution. The first, second, and 
third peaks are at - 560, 890, and 1120 MeV/c 2 , re- 
spectively. The first and second peaks have very little 
background. Fits to the first two peaks using Gaus- 
sians for the signal shapes result in peak positions of 
(561.7 ± 0.1)MeV/c 2 and (893.4 ± 0.2) MeV/c 2 . These 
values are compatible with the (ra-f(2S) ~ m T(is)) an( i 
(m X ( 3S ) - m T (is) ) PDG values, respectively. We con- 
clude that the first and second peaks are due to the de- 
cays T(2S) -> T(lS')7r+7r- and T(3S) -> T(lSV + 7i-, 
where the T(2S, 3S) are produced mainly in the radiative 
return processes e + e~ — > T(2S,3S)j. 

In contrast to the first two peaks, the third peak has 
modest background. The position of the peak is derived 
from a fit to the AM distribution in the third peak re- 
gion (see Fig. \2§ using a Gaussian for the signal and a 
third order polynomial for the background. The result is 
AM = (1119.2 ± 0.4)MeV/c 2 , which is in good agree- 
ment with the mass difference (mx(4S) — m T(is)) from the 
PDG. The Gaussian width is a = (5.7 ± 1.0) MeV/c 2 , 
which is consistent with the estimated AM resolution. 
The signal yield in the interval 1110McV/c 2 < AM < 



nal yield fixed at zero and its best fit value, respec- 
tively. This peak is identified as a signal for the decay 
T(45) -> T(lS)7r+7r- with a subsequent T(15) -► 
transition. 

To verify this interpretation, we study the resonance 
properties in more detail. First, the AM distribution 
for the (/i + /i~7r + 7r~ X) event sample in the M M +^- > 
8 GeV/c 2 mass region is considered. Here we use a looser 
requirement on M^+^- to study the background in a 
wider region. The distribution of AM vs M„+„- for 
the on-resonancc (off-resonance, ^/s = 10.52 GeV, inte- 
grated luminosity J Cdt = 49.4 fb _1 ) sample is shown 
in Fig. [3^(b). Comparing the two distributions, we see 
that the behaviour for the on-resonance and off-resonance 
sample are similar except that in the off-resonance data, 
the third peak is absent. The "slanted band" in this plot 
is due to the non-resonant reaction e + e~ — > fi + fi~n + Tt~ . 
Using the off-resonance sample, the AM distribution 

m T (is)| < 



where M^+^~ is restricted to |M M + A 



60MeV/c 2 has only two peaks, corresponding to T(25) 
and T(35) decays. For the on-resonance sample when the 
data size is scaled to the off-resonance sample, the total 
number of events and the number of background events 
in the interval 1110MeV/c 2 < AM < 1135MeV/c 2 
is = (7.3 ± 0.9) and A^ s g = (2.2 ± 0.5), re- 

spectively. If we consider a AM interval shifted by 
— 60MeV/c 2 , to take the lower yfs into account, we find 



1135MeV/c 2 determined from the fit is A^ cv = 43.9 ± ^tot = (3.0 ± 1.7) events in the off-resonance sample. To 
7.9, with a statistical significance of 8.0 a which cor- compare the scaled total (N t r ° t s ) and background (iV^ s g ) 
responds to — 21n(£ /£max) = 72.3 with 3 degrees of 
freedom (mass, width, and yield). Here Cq and £ max 
are the likelihood values returned by the fit with sig- 



yields as models for the off- resonant yield N°^ t , we use 



the so-called likelihood x 
N hkg> and X 2 = 3.26 for N] 



res . 
tot ' 



We find x 2 = 0.26 for 
the background interpre- 
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FIG. 



3: The AMvsM p+A 



scatter plot for the on-resonance sample (a) and off-resonance sample (b) in the M M + M - > 8 GeV/c 2 
mass region. The rectangle in (a) indicates where T(45) — > T(lS)7r + 7r _ transitions would cluster; the qq background in this 
region should scale from that in the rectangle in (b). For clarity, the rectangle in (a) is drawn larger than the region used 
to determine the signal yield (which is determined by optimizing the expected signal yield relative to the square root of the 
expected background). 



tation of the events in the third peak region is therefore 
favoured, although the discrimination is relatively weak 
(1.7a). This suggests that the events in the third peak 
region show no sign of a similar enhancement. 

Additional information can be obtained from the 
study of the 7r + 7r~ system. The efficiency-corrected 
distribution of the invariant mass of the 7r + 7r~ system 
(M v +„-) is shown in Fig. [4] for events in the third 
peak region. The background (see Fig. [2]) is not sub- 
tracted from the M„+^- distribution. The efficiency 
is calculated by a Monte Carlo simulation. The Evt 



Gen event generator [101 ] with a matrix element 11 1 
taking into account particle spins is used to produce 
T(45) — * T(lS)7r + 7r _ — > /i + /i~7r + 7r~ events that are 
passed through the detector simulation and reconstruc- 
tion programs. 

As shown in Fig. 2J the M n + n - distribution for events 
from the third peak region in the AM distribution shows 
an enhancement at high masses. In contrast, the distri- 
bution for background events, which is taken from the 
AM sideband (1004MeV/c 2 < AM < 1110MeV/c 2 , 
1135MeV/c 2 < AM < 1210MeV/c 2 ) and normalized 
to the background events underneath the peak, is more 
uniform and shows no sign of similar enhancement. This 



difference in the behaviour of the M„ 



distribution is 



an additional argument in favour of a resonance interpre- 
tation of the third peak. 

The M^^- distribution can be described by the shape 



predicted by the Brown-Cahn model [ll| (see Fig. 4), 
where the hadronic transition between heavy quarkonia is 
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FIG. 4: The tv + ty~ invariant mass distributions for events 
from the third peak region in the AM distribution. The 
shaded histogram is for background events estimated from 
the AM sideband. The solid line shows the AM distribution 
predicted by the Brown-Cahn model [ill ]. 



considered as a two-step process: the emission of gluons 
from heavy quarks and subsequent conversion of these 
gluons into light hadrons. 

The branching fraction for the T(45) -► T(lS , )vr+7r- 
decay is determined from B(T(4S) — > T(lS')7r + 7r _ ) = 
^Vcv/(^Vt(4S) • £ ■ B{T(\S) -> fi+fJ,-)). The total num- 
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ber of T(45) in the data sample is Nxus) = (464 ± 
6) x 10 6 , and the nominal branching fraction B(T(1S) — ► 
= (2.48 ± 0.05)% 4.]. The efficiency obtained 
from the Monte Carlo sample is s = (2.14 ± 0.06)%. 
We apply a correction of about 8% to E sum /^/s, one of 
the variables used for selecting hadronic events. The ef- 
ficiency is sensitive to this variable; the correction is de- 
signed so that the data and MC match in this variable. 

The systematic error of the reconstruction efficiency 
due to this correction is 8%. The systematic uncertainty 
in the reconstruction efficiency due to lack of knowledge 
of the T(AS) — > Y(lS)ir + ir~ — > /i + /.i^tt + tt~ decay ma- 
trix element is estimated by varying the parameterization 
of the M n + V - distribution. Requiring a reasonable fit to 
the M„.+„.~- distribution we find a 3% variation in the ef- 
ficiency. Other systematic uncertainties come from the 
choice of the fit range and the background shape (4.5%) 
in the AM distribution, choice of the signal range (4%), 
choice of the T(15) mass range (6%), and from the track- 
ing efficiency (4%). The total systematic uncertainty is 
obtained by adding these contributions in quadrature; 
the result is 13%. 

The measured product branching fraction is 

B(T(AS) -» T(lS>+7r-) x B(T(LS) -» n+n~) 
= (4.4 ± 0.8(stat.) ± 0.6(sys.)) x 10~ 6 . 

The branching fraction is 

B(T(45) T(lS)7r+7r-) 

= (1.8 ± 0.3(stat.) ± 0.2(sys.)) x 10~ 4 . 

We also extract the partial decay width for 
T(AS) — > T(lS)n + Tr~ transition using the world- 
average value of the total width [l| and obtain 

r(T(45) -> T(lS , )vr+7r-) 

= (3.7 ± 0.6(stat.) ± 0.7(sys.)) keV. 

The measured values of B(T(45) -> T(15)7r+7r-) and 
r(T(45) — *■ T(1S)tt + tv~) are about twice larger than 
BaBar's results Q. 

To summarize, a study of transitions between T states 
with the emission of charged pions has been performed at 
Belle. The mass difference distribution {M ll + li - 7T + 7I - — 
M^+p-) from the fi + 7r + X event sample for M M + M - 



within the T(IS') mass region has two peaks from 
T(25, 3S) — > T(15)7r + 7r _ decays with no background. 
A third peak at AM = (1119.2 ± 0.4) MeV/c 2 is inter- 
preted as a signal for the decay T(45) — > T(lS)ir + ir~ 
with a subsequent T(15) — * /i + yT transition. This final 
state is the first example of a non-BB decay mode of the 
T(45). The branching fraction B(T{AS) -> T(15)7r+7r-) 
and the partial decay width T(T(45) — > T(l,5')7r + 7r _ ) are 
measured. The M^^- distribution can be described by 
the shape predicted by the Brown-Cahn model fill ]. 
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